BB mixing with the bulk fields in the Randall- Sundrum model 



Jong-Phil Lee* 

Department of Physics and IPAP, Yonsei University, Seoul, 120-749, Korea 



Abstract 



We calculate the BB mixing in the Randall- Sundrum bulk model. In this 
model, all the Standard Model fields except the Higgs can reside in the bulk. 
Two suggestive models of "mixed" and "relaxed" scenarios are considered. 
We find that the enhancement of the loop function is 0.51% for the "relaxed" 
and 1.07% for the "mixed" scenario when the first 4th KK modes are included, 
for a bulk fermion mass parameter i' = —0.3. 
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I. INTRODUCTION 



The idea of extra dimension provides a very elegant explanation for the gauge hierar- 
chy problem of the standard model (SM). Arkani-Hamed, Dimopoulos, and Dvali (ADD) 
suggested that there exist n fiat extra dimensions with factorizable geometry [1]. In this 
model, the enormous Planck scale Mpi is induced by the largeness of the extra dimensional 
volume Vn, via Mj,i = Mo"+V„. New fundamental scale Mq in 4 + n dimension can be set 
down around 1 TeV, which resolves the hierarchy problem. An alternative to ADD was soon 
proposed by Randall and Sundrum (RS) [2] where the gauge hierarchy is explained by an 
exponential warp factor from a 5-dimensional non-factorizable geometry. In this model, two 
branes are embedded at the boundaries of the AdSs slice with a single S1/Z2 orbifold extra 
dimension. 

In the original model of ADD or RS, only the gravity can propagate in the bulk. Effects 
of the bulk graviton appear in 4D world as an infinite tower of Kaluza-Klein (KK) modes. 
Their phenomenological signatures at colliders are widely studied. 

A natural extension at the next step is to put the SM fields into the bulk. In the universal 
extra dimension (UED) model by Appelquist, Cheng, and Dobrcscu (ACD), all the SM fields 
are allowed to live in the extra dimensions [3] . In the RS model, Goldberger and Wise tried to 
put the scalar fields into the bulk, which has been developed to the bulk-stabilizing modulus 
field, or radion [4]. They provided an open possibility that the SM matter fields might 
reside in the extra dimension. Later the gauge bosons are placed in the bulk in [5]. But the 
electroweak precision data constrain the first KK state of the gauge boson so strongly that 
the typical scale on the TeV brane goes up to about 100 TeV. Apparently, this is not a good 
solution to the gauge hierarchy problem. The problem was alleviated by putting fermion 
fields in the bulk while confining the Higgs to the TeV brane [6] . An attempt of constructing 
the bulk fermion fields was already done in [7] to explain the neutrino masses and mixings. 

Equipped with the RS-bulk SM fields, lots of works have been done to refine the model 
and accommodate it to the existing data. Here the fermion mass parameter u = m^/k 
plays an important role, where is the fermion bulk mass and TZ^ = —20k^ is the 5- 
dimensional curvature [8,9]. When placing the fermion fields in the bulk, there might occur 
large contributions to the flavor changing neutral current (FCNC) or the SM p parameter. 
Simple assumptions of a universal bulk fermion mass and the minimal flavor violation where 
the CKM paradigm governs the flavor mixing can be a solution to this potential problems. 

However, when the Yukawa interaction is taken into account, non-negligible mixing in the 
top sector can give a large contribution to the p parameter. Even worse is that the shifted 
mass spectrum in the top quark KK modes does not guarantee the Glashow-Iliopoulos- 
Maiani (GIM) cancellation, which will cause a disastrous FCNC. Toward a more realistic 
model, Hewett-Petriello-Rizzo (HPR) proposed the so called "mixed scenario" where the 
first two generations of fermions are placed in the bulk while the third is localized on the 
TeV brane [10]. While the HPR model reproduces the quark mass hierarchies nic/mt and 
ms/mf, without spoiling the p parameter constraints, it seems unnatural to confine only one 
family into the wall; there still remains a potential danger of FCNC. 

As an alternative, a "relaxed" model was suggested by Kim-Kim-Song (KKS) where the 
assumption of universal bulk fermion mass is slightly modified, retaining all the fermions in 
the bulk [11]. In this approach, the S'[/(2)-singlet bottom quark field has a different bulk 
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mass m'^. They showed that introducing another parameter can be accommodated well to 
the electroweak precision data of Ap = p — PsM and b S7. 

In this paper, we analyze BB mixing based on the KKS. With the advent of the B- 
factory era, we are now entering the age of precision tests of flavor physics. Current world 
average of AM^, which parameterizes the BB mixing, is AMb — 0.502 ± 0.006 ps~^ [12] 
where the experimental error is very small. Relevant box diagrams are already studied in 
ACD-UED model, resulting in a 17% enhancement of the loop function [13,14]. Present 
work will be a good comparison to the ACD-UED result. It is also interesting to compare 
the results from KKS and HPR each other. If the fermion mass parameter u — —0.3 is 
chosen (which is phenomenologically viable in the literature), especially, mixing in the top 
sector is rather small; the GIM cancellation is incomplete but we expect some remnant. And 
the diagonalization of the mass matrix can be done perturbatively. Naive thoughts lead to 
the idea that HPR will produce larger effects than KKS, yielding to what amount of the 
excess of FCNC from the HPR would appear. 

The paper is organized as follows. The RS-bulk model is reviewed in the next Section. 
We simply omit the gravity in the set up. In Sec. Ill, mixing in the top sector is considered. 
We follow the description of KKS [11], and do not consider possible gauge boson mixing 
for simplicity. The box diagrams for the BB mixing are calculated in Sec. IV. Section V 
contains the results and discussions. A summary is given in Sec. VI. 



II. SETUP OF THE MODEL 

In the RS model, one spatial dimension is compactified on a S1/Z2 orbifold of radius Tc 
with a nonfactorizable geometry 

where the four-dimensional metric tensor is rj^y = diag(l, — 1, — 1, — 1), and cr(0) = 
krc\(f>\, < [01 < TT. The model parameter k is related to the 5-dimensional curvature 
7^5 = -20A;2. 

In the original RS model, only the graviton can propagate through the bulk with its KK 
modes. We do not consider the graviton KK modes here because they are irrelevant for BB 
mixing. We follow the model setup of KKS [11]. 

Nonabelian gauge fields A'^{x, (p) can reside in the bulk via the 5D action 

Sa^-IJ , (2) 

where = 9m - OnA^i - g^e^-^'^A^j^^A^ {a,b,c = 1,2,3). Choosing the gauge of 

A4{x, (/)) — and assuming the KK expansion of Af^ to be 

A«(x,0)=5]Af)(x)^^, (3) 

n=0 V'^c 

we have a 4-dimensional effective action of massive KK gauge bosons as 
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/°° r 1 1 
n=0 L 



(4) 



if the extra dimensional component Xa' (^) the Bessel function 



(n) e' 



^(0) r 



N 



(n) 



in). 



Here 



k 



(5) 



(6) 



where M^"^ is the mass of the nth KK mode of the gauge boson. The continuity of dxA^ /d(j) 
ai (f) — and = ±7r determines the mass spectrum and coefficient , 

,(n)_ Ji(My/fc) + (Mi"Vfc)J{(Mi"Vfc) 



(at = 0), 



Ji(4")) + J((4"^) + + =0 (at = ±7r) , 

where x^^^ = z^^\(l) = tt) = (M^^ /k)e^'^''^ . The normahzation constant 



(7a) 
(7b) 



J,(4")) + a?V,(.J)) 



4"'(^=0) 



is obtained by the orthonormahty condition 



/7i 
■7 



(8) 



(9) 



Now consider the bulk fermion with arbitrary Dirac bulk mass. A Dirac fermion field ^ 
with bulk mass is described by the 5-dimensional action [7] 



- m^sgn(0)\l'\l' > , 



(10) 



where Da is the covariant derivative, 7— = (t'^j^Ts), and the inverse vielbein E-^ — 
diag(e'^, e'^, e'^, e'^, 1/rc). The underlined uppercase Roman indices describe objects in the 
tangent frame. Possible spin connection cubca is omitted because it vanishes when the 
Hermitian conjugate is included. 
Integration by parts leads to 



1 

'2^ 



(11) 



where the periodic boundary conditions of '^L,R{x,'!r) — '^l,r{x, — tt) are imposed. Expand- 
ing ^ as 

*,,«(x,0) = ^V^a^/S('^)' (12) 

n=0 V^c 

and requiring that 

r #e-('^)ii-)*(0)/f (0) = r #e'^(^)/ir^*(0)/i,"V) = 5- , (13) 



' r. 



(14) 



(15) 



we have the effective action for the massive Dirac fermions 

Sf^Y. [ d^xUp^''\x)ip iP^''\x) - M)"^^(")(x)^(")(a:) 

n=0 

The Za-symmetry of Sp in (10) requires that = '^l'^r + '^r'^l be Za-odd. This 
requirement is satisfied if we impose the opposite Z2-parity on /|"^. We fix f^j^^ to be 
Z2-even(odd). Introducing u = m^/k of order 1, the solutions are, for n 7^ 0, 



/(")^r(")(0) 



Ji/2-.(^^'^))+/5j^^l^i/2-.(^("^) 

Jl/2+.(^('^^)+/5^^n/2+.(^("^)' 



where zf\^) = {Mf^ /k)e<'^^ and for n = 0, 



r(") 



(16a) 

(16b) 
(16c) 



(17) 



Mass spectrum and the coefficients jS^}- as well as the normalization constants N^)- are 
determined by the boundary conditions: 



d 



-(") 



</>=0,7r 



(18) 



Explicitly, 



J_(i/2+.)(m}"VA:) 
yi(i/2+.)(M}"V^) 
from the boundary conditions at = 0, and 



^(n) 



(19) 



5 



(20) 



at = TT, where x^^^ = z'^\ 



tt). Similar expressions can be given for the t^"^(0) sector, 



and the left- and right-handed excitation masses are degenerate to m|"'* for a given n. The 
normalization constants are 



X,r 



(0=0) 



(21) 



Some difficulties arise when the fcrmions are included in the bulk because the fermion 
field contents must be doubled. In our setup of (11), a fermion which belongs to a specific 
representation of a gauge group should be vector-like, possessing both left- and right-handed 
chiralities. For each generation, we have an S'[/(2)-doublet Q — {qwQ.dY and two SU{2)- 
singlets u and d where both chiralities are allowed for all of them. Exphcitly, 



Q{x, 0) = + = ^ 

n 



u 



Qt\x)x^-\<P) + Q^^\x)t^-H<P) 
^''\x)T^^\ct>)-rut\x)x^^\ct>) 



(22a) 
(22b) 



where we have assigned x*^"^ to the left-handed SU{2) doublet and the right-handed SU{2) 
singlet in order to accommodate the SM fermion fields. The SU (2) singlet down part d{x, 0) 
has the same KK decomposition as u{x, 0). 

With these bulk-SM fields, the charged current interactions look like 



Sffw = J d'xV^e''^[qu\f^qd + h. 



c. 



/ 



v ^ /=0 '-n,m=0 



n,m=l 



+ h.c. 



(23) 



where g — g^s/v^Trr^ and the coupling constant is 



nml 



e'^x(")(0)x(™H0)x?(0) ■ 



(24) 



III. MIXING IN THE TOP SECTOR 

For the bulk fermions, there are two sources of physical fermion mass spectrum. One 
is the KK mass which is proportional to ~ QulQuR or ~ ulUr. The other is the Yukawa 
interaction which relates the SU{2) doublet with the singlet, ~ QulUr, as in the SM. This 
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difference results in mixing among tlie fermion KK modes. If tlie quark masses were suffi- 
ciently smaller than the KK mass scale, this mixing effects might be negligible. In this case, 
all the fermion KK mass spectrum would be almost degenerate assuming that the model is 
minimal with a single parameter for the universal bulk fermion mass. The exact degen- 
eracy between the KK masses for the T3 = ±1/2 fermions is responsible for the vanishing 
contribution to Ap. Furthermore, wc can expect the GIM cancellation at each level of KK 
modes because all the up- (and down-) type quarks arc degenerate. 

The presence of top quark which is very heavy makes the mixing non-negligible. Con- 
sequently, degenerate top and bottom KK masses are shifted. They can cause a large Ap, 
because the quantum correction is proportional to the squared mass difference between up- 
and down-type quarks. Adding higher KK modes makes the problem worse. 

As mentioned in the Introduction, HRP proposed the "mixed" scenario where the 3rd- 
generation fermions are confined to the TeV brane while the other two can propagate in the 
bulk. Instead, we adopt the KKS model of [11] where the 5'C/(2)-singlet bottom quark field 
has a different bulk fermion mass m^. The model is recapitulated in the 5D action for the 
3rd-generation quarks 



-sgn(0)(m^{(5Q + tt} + m'J)h) 



(25) 



Expanding the fermion fields with the mode functions and integrating over yield the bulk 
fermion KK masses as 



00 



n=l 



+ h.c. , 



(26) 



where kEw = ke = ke''^^"'^, and u' = m'^/k. 

Now consider the Yukawa interactions associating the Higgs boson. Ordinary Higgs 
mechanism works here so that the SM particles acquire masses. It has been argued, however, 
that if there were bulk Higgs fields, the hierarchy problem remains unsolved [6] or the 
observed W and Z mass relations cannot be reproduced [9] . Wc assume here that the Higgs 
field is confined to the TeV brane. Then the 5D action for Yukawa interactions is 



^Q{x, 0) • H{x)h{x, 0) + ^e'^'Qix, 0)„ • H{x),t{x, 0) + h.c. 
k k 



5{4> - tt) , 
(27) 

where A5'* is the 5D Yukawa couphng. After the spontaneous symmetry breaking, the Higgs 
field shifts ^ v^-\- H'^, and the 4D effective Lagrangian becomes 



£e« = ^(^;?+xi.ii^+--o(4°^+xi4^^+---) 



(28) 



where = Xf{l + 2v)/2{l - e'+^''), v = ev^, Xn = X^^'K^, ^)/x^''K^, and x'n are with 
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Prom Eqs. (26) and (28), the physical mass term for the top sector is 



f(0) 



V '■■ } 



(29) 



The bottom sector shows very similar mass terms. Let the number of the KK states be n^o 
which is, in principle, infinite. Then the i^n^ + 1) x (2noo + 1) matrix }At has the form of 



M% 



(30) 



where the (^oo + 1) x (no© + 1) matrix My is from Yukawa mass term and the (^oo + 1) x rio 
matrix A4kk from KK masses. They are given by 



/ 1 Xl X2 

xi xl X1X2 ■■■ 

X: 

V ^ 



X2 X2X1 X2 



\ 



(31a) 



/O 



X 



(1) 



a; 





(2) 
/ 





V: : : 



KK 1 



(31b) 



where rritfl = Xtv/\/2. 

Now the physical mass eigenstates u'^^^ are obtained by diagonalizing (30) through an 
orthogonal matrix J\f: 



/ '(0) \ 



/(2) 



.2 



u 



(1) 



(32) 



V : / 



For example, q^^}^ can be written in terms of KK mass eigenstates 



2no 



j=0 



'(i) 



(33) 
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Hereafter, the underlined index runs from zero to 2noo- In terms of mass eigenstates, the 
Wtb-vertex in 4D becomes 

^ = ^1^*" ( E ] ypYt^WlP + h.c. (34) 



\m=0 



For hght quarks (m^o = 0), A4q can be analytically diagonalized to 
diag{mP , Mi°\ , ■ ■ •) where 



ml 



f=0, Ml°) = = W . (35) 



The corresponding orthonormal matrix M"^^^ is 



•'^(0,0) — ' •'^{0,n) — -'^(n,0) " ^ ' 

a/2 ■ 



On the other hand, the diagonalization of M.t of the top sector is nontrivial since rritfi 
is heavy. However, an approximate calculation is possible unless Xn'^ much larger than 
unity, with a small parameter rritfl/kEw- K is shown in [11] that a perturbative diagonal- 
ization can be done for v > —0.3. Up to the leading order of the expanding parameter 
5 = rritfi/kEw, the mass eigenvalues are 



mg = mo , Mn = kEw (^xf^ + 

^2 



Mn^+n = kEw [xf^ - , (37) 

and the orthonormal matrix J\f are 

■A/'fezr.)^<L)+<L);|' (38) 



where 



2 



'0) ~ ^(") ' 

(n,n) •' {noo+n,n) 4-^(") ' 



(n) , . 

<l,-<i«,™,--=^^fe?^ (forn^m). (39) 



9 



IV. BB MIXING AND THE BOX DIAGRAMS 



The BB mixing is parametrized by the physical mass difference 

AMb = ^\mn,s{AB = 2)\B')\ 
Mb 

= ^^iVtlVtdl'BBflMsVBSoixt) . (40) 

where Mb (Mw) is the B (W) mass, r]B the QCD corrections, and Bb and Jb are the bag 
parameter and decay constant, respectively. The loop function SQ{xt) which recapitulates 
the box diagrams is 

ixt-llxt + x^ Sx^lnxt , . 

- 4{i-x,r 2(1^ ' 

in the SM, with xt = {nit/MwY- 

In the RS-bulk model, bulk fermions and gauge bosons contribute to the box diagrams, 
as shown in Fig. 1. Since the external fermions are the zero modes, relevant vertex factor is 
CfJ(^ multiplied by the mass-diagonalizing orthonormal matrix J\f. 

Introducing the effective couplings 

T'^ = (at . C^('i')t^y\ (42a) 
Q'^ = ■ C^^'^^'^^y , (42b) 
the loop function has the form of 



S(m,n;i,j) = M^J^ 



1 + "^'^ ^^'^ „ Q^"^Q^"Q^"Q^"'/(W^('")W^(")g(*)g(^)) 

4M^") Mi") / 

(m^(,)m^(,, + m^ajm^,,, \ , ^ ,,,,,, 

4Mj")'Mi")' ) 

+ I + 1 1 - m^«m,'o)r"'r"T^"r^"*7(iy("V^"')^^')^^')) 



+ ' 



Here we separate the longitudinal part of W^'^\ denoted by for the calculational con- 
venience, and 
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Jo 

/(V^(n)0(m)^«^O)-) ^ f ^^^y 

Jo 

where 

mli),u) = ^rnl^ + (1 - ^)m]u) , <n = yM"^' + (1 - ■ (45) 

V. RESULTS AND DISCUSSIONS 

For the numerical results, we used kvc — 11.5, ksw = 1 TeV, and u = —0.3. The 
choice of z/ = —0.3 allows us to use Eqs. (38) and (39); for u < —0.4, only the numerical 
diagonalization of the mass matrix is reliable. One question may arise whether this choice 
would spoil the Ap accommodation. This, however, would not be the case because we 
have another adjustable parameter u' for the 6-sector in KKS model. Even in the case of 
(i/, I/') = (—0.3, —0.6) the mass differences between top and bottom KK modes are much 
smaller than 1 TeV [15]. We also give the results for u = —0.2, — 0.1 for comparisons. 

In this analysis, we include the KK modes only up to n = 4 gauge boson (which corre- 
sponds to n = 2 fermion because bulk fermion has double field contents) for simplicity. As 
a comparison, the HPR contribution to the box diagrams is also given. Numerical results 
are summarized in Table I and Figs. 2, 3. Main results of the analysis are as follows. 

First, we have the loop function enhancement 0.51% (1.07%) for u = —0.3 in KKS 
(HPR) up to n = 4, as shown in Table 1 and Fig. 2. The numbers can be compared to the 
ACD-UED result where the increase amounts to 17 % [14]. Increase of the loop function 
imphes the smaller value of \Vtd\ , because AMb — 0.502 ± 0.006 ps~^ ~ l^dP ^ S{m, n; i, j) 
is now well measured [12]. Figure 3 shows our results for \Vtd\- We used a/ Bsfs — 235 
MeV, r]B — 0.55. The QCD correction factor riB has a small error ±0.01 compared to the 
hadronic uncertainty [14,16]: 

/^/b = (235il?) MeV . (46) 

Since the BgBs mixing shares the same loop function with BB mixing, increase in 
Sim, n; predicts the equal amount of larger AMb^- Current bound of it is AMb^ > 
15 ps-^ [12]. 

In order to strongly constrain the model parameters, or even rule out, measurements 
of sin 2/3, e.g., would not be crucial. Rather, determining the CKM angle 7 as well as 
reducing the hadronic uncertainty significantly will severely restrict the parameter space. 
For example, in the ACD-UED model, Ref. [13] argued that larger value of a/ BbJb above 
222 MeV with one-third reduced error will rule out the the UED model or push up the 
compactification scale up to multi-TeV region. An analytic study is much easier in the ACD- 
UED model compared to the RS-type one because the mass spectrum and mode functions 
have simpler structure in the UED model [3]. In the RS-bulk models, on the other hand, 
the n-dependence of the mass spectrums, e.g., is not explicit. 



nr 



2rrr 



Ml 



mn 



In 



m 



m 



(<)oO); 



q(i)qU) 
mn 



, (44a) 



(M, 



2 

mn 



+ 



m 



qy-'q' 



mn 



In 



m 



m 



qy-'q 



M2 
mn 



, (44b) 
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Second, the HPR gives larger values than KKS for each n. The reason is as follows. If 
all the fermions were in the bulk and no mixing happened, there would be GIM cancellation 
at every n-th excitations. In case of KKS, the cancellation is incomplete because there is 
a small mixing (when u > —0.3) in the top sector. In HPR, however, the third generation 
resides in the 3-brane; only the first and second generations have their KK excitations in the 
bulk. There is no compensating contribution from the third generation KK modes, which 
results in larger loop function. This feature was already pointed out in [11]. 

The GIM cancellation is advantageous in viewpoint of the convergence. There is, of 
course, no way to warrant the finiteness of adding up infinite KK tower by checking just 
a few excitations. In the ACD-UED, it is argued that the GIM mechanism improves the 
convergence of the loop function in BB mixing [14]. And it is quite encouraging that the 
RS-bulk SM effects of KKS on 6 — > 57 are satisfactorily converging as we add higher KK 
modes [11]. 

Third, the i/-dependence is shown in Fig. 2 and Table I. The KK mode contributions 
to the box diagrams get larger as u grows. The reason is that both the coupling and the 
bulk fermion mass become larger when u goes from —3 to —1. This pattern is already 
studied in previous works; see Fig. 1 of [9] and [10]. We have, for example, the bulk fermion 
masses xf^EW = 2.71, 2.85, 3.00 (TeV) for u = -0.3, - 0.2, - 0.1, respectively. It 
is already suggested that the plausible range ofi/ is —0.8 < u < —0.3 in the literature 
[9,10]. If > —0.3, constraints from the electroweak precision data make the first KK mode 
of gauge boson heavier. In our case, since the loop function grows abruptly for large u, 
present analysis of BB mixing can provide another hint for the upper bound of v when the 
uncertainties of the parameters for AMb are reduced, just as in [13]. 

One thing to be noticed is that BB mixing involves only the up-type quark's mass 
parameter u. The main motivations of KKS are i) all the fermions can reside in the bulk, ii) 
different mass parameters for the top- and bottom-sector can accommodate the electroweak 
precision data. In this point of view, DD is a good testing ground for the bottom sector. 
A preliminary result of our approach to DD mixing is that the contribution of the first few 
KK modes is extremely small. This is mainly due to the smallness of the b quark mass. 
Consequently, mixing with the KK masses is quite weak and the GIM cancellation is more 
complete. Unfortunately, current experimental results are rather poor [17]. More reliable 
check on the model can be implemented by simultaneous analysis on BB and DD with 
improved measurements, as well as on the electroweak precision tests. 

VI. SUMMARY 

In this paper, we have calculated the BB mixing in the RS-bulk model. Contributions 
from the newly introduced KK modes are positive to increase the loop function. We found 

that KKS predicts smaller contribution than HPR due to the remnant of the GIM cancella- 
tion. More precise determination of the CKM unitarity triangle will test the validity of the 
model. Also, further developments of the present work, such as including much more KK 
modes or simultaneous application to various observables, remain challenging. 
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FIGURE CAPTIONS 



Fig. 1 

Bulk gauge bosons and bulk fermions contribute to the box diagrams. 
Fig. 2 

Enhancements of the loop function ^S{m,n;i,j) for (a) KKS and (b) HPR for various u. 
The X-axis is the number of KK modes included. 

Fig. 3 

Results for \Vtd\ from the loop functions for KKS and HPR. Each hues corresponds to 
u — —0.1, —0.2, —0.3 from bottom to top, respectively. The top line is the SM value. 
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TABLE CAPTIONS 

Table 1 

Enhancements of the loop fuction ^ S{m, n; i, j) in % for KKS and HPR for different values 
of v. 
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TABLES 



TABLE L 





u = -0.3 


u = -0.2 


u = -0.1 


KKS 


0.51 


2.06 


26.8 


HPR 


L07 


3.98 


50.5 
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